Nanotube-Metal Junctions: 2- and 3- Terminal Electrical Transport 
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We address the quality of electrical contact between carbon nanotubes and metallic electrodes by performing 
first-principles calculations for the electron transmission through ideal 2- and 3-terminal junctions, thus reveal- 
ing the physical limit of tube-metal conduction. The structural model constructed involves surrounding the tube 
by the metal atoms of the electrode as in most experiments; we consider metallic (5,5) and ;i-doped semicon- 
ducting (10,0) tubes surrounded by Au or Pd. In the case of metallic tubes, the contact conductance is shown 
to approach the ideal 4e 2 /h in the limit of large contact area. For three-terminals, the division of flux among 
the different transmission channels depends strongly on the metal material. A Pd electrode has nearly perfect 
tube-electrode transmission and therefore turns off the straight transport along the tube. Our results are in good 
agreement with some recent experimental reports and clarify a fundamental discrepancy between theory and 
experiment. 
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Carbon nanotubes (CNTs) have received extensive experi- 
mental attention for more than a decade 12 - 3 , and are consid- 
ered a possible basis for nanoelectronic technology indepen- 
dent of silicon. A major issue is the quality of CNT/metal 
contacts: obtaining the minimum contact resistance is critical 
to access the intrinsic electric properties of CNTs. Despite 
extensive experimental effort to improve the contact trans- 
parency and reveal the relevant factors behind it - metal mate- 
rial, contact structure, and type of tube, for instance - a clear 
picture is still not available. 

On the theoretical side, it is highly desirable to be able 
to simulate from the first principles the electron transport 
through CNT/metal junctions and thus to improve our under- 
standing of this important issue. So far, first-principles studies 
of the contact transparency between metallic CNTs and metals 
have been carried out [(3,3), (4,4), or (5,5) tubes with Al, Au, 
or Ti, for instance The results are, however, quite scat- 
tered, and agreement between theory and experiment has not 
yet been achieved: For Al electrodes, an equilibrium conduc- 
tance of ~ 1 Go (= 2e 2 /h, the conductance quanta) was found 
by some calculations 4 ' 6 while ~ 2 Go was found by another 5 . 
For Ti electrodes, 1.7 Go was obtained by one calculation 7 
and ~ 1.2 Go by another 6 , while experiment 8 found ~ 2 Go- 
For Au electrodes, a recent calculation 6 showed ~ 1 Go while 
a value of 1.5 Go was found experimentally 9 . On the other 
hand, a model calculation 10 using a jellium model for the elec- 
trode even showed that the ideal CNT/metal conductance will 
not be larger than 1 Go- With regard to the metal used, it was 
found theoretically that Ti is better than either Au or Al for 
contact transparency 6 . On the other hand, experimentally, Pd, 
for which no theoretical calculation is available, was found 
superior to Till. 

It has been unclear how to explain these discrepancies. One 
possible reason, as is supported by this work, is the differ- 
ence in contact structure between theory and experiment. In 
all these calculations, simplified contact models were adopted 
due to computational cost or methodology: On the carbon 
side, the tube-metal connection is made by either straight u- 
bonds from the tube end or ir -bonds from the side of the tube, 
while the metal is modeled by a thin nanowire or small clus- 
ter. However, in most experimental situations, a tube is sur- 



rounded by metal atomsi^iiiiii 5 ^! 6 *!!. As the contact structure 
and quality changes from case to case, experimental results 
are also scattered. In such a situation, the key contribution 
to be made by theory is to reveal the physical limit to which 
experimental measurement may approach by improving the 
contact quality. 

Furthermore, compared to metallic CNTs, semiconducting 
CNTs are much more important for potential electronic appli- 
cations. This naturally raises the issue of contact transparency 
between doped semiconducting tubes and metal junctions. In 
addition, electron transport through multi-terminal structures 
is a key property in moving toward applications. Neither of 
these fundamental issues has been tackled previously, to the 
best of our knowledge, using first-principles calculations. 



In this paper, we first construct a betteji^iiiiiii^ji&iii struc- 
tural model of an ideal contact, as shown in Fig. \J]for the 
device region, in which an infinitely long tube is surrounded 
by metal atoms of an electrode. We use this structural model, 
first, to investigate the two-terminal tube-metal transparency, 
Pterin, (from terminal 1 to 3), obtaining the physical limit of 
tube-metal conductance for an ideal contact. Second, we then 
study three-terminal electron transport by including explicitly 
the self-energies of all the leads shown in Fig. ^ Specifi- 
cally, we investigate the division of current among transmis- 
sion straight-through the junction, probability to turn the cor- 
ner, and reflection. In both 2- and 3- terminal cases, we con- 
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FIG. 1: Side and top views of the device region of a (5,5)-tube-Au 
junction, where the Au electrode consists of 7 atomic layers (7L). 
The length of the tube in the device region is 32. 6A and the dimen- 
sions in the lateral directions are 20.4 and 22.44A, respectively. The 
three terminals are denoted by ti, ta, and t3. 
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FIG. 2: Two-terminal transmission T2tcrm. (E) (tube to metal) for (a) (5,5)-Au and (b) (5,5)-Pd with different widths of the metal electrode, as 
indicated in the legends. Also plotted is the transmission for a pure (5,5) tube, as a comparison. 



sider the dependence on the contact quality (simulated by the 
width of the metal electrode), the metal material, and the type 
of tube (metallic or semiconducting). 

The structural model (Fig. ^ can be regarded as an ideal 
CNT- metal junction because of three features: (1) The tube is 
perfect and infinitely long, thus eliminating the quantization 
in the transport direction which bedevils calculations 6 using 
a short CNT segment. (2) The tube is surrounded by metal 
atoms, as is usual in experimentsiSililiil^iiki!. (3) The large 
active region around the countact ensures that the tube-metal 
interaction is fully included. 

Contact quality will depend experimentally on the number 
of good carbon-metal connections, which is essentially deter- 
mined by whether the metal wets the CNT surface. If the 
metal (like Au) does not wet the CNT surface, there will be 
few connections even though the electrode is large. In our 
models, because of the ideal structure, the number of good 
carbon-metal connections is substantial even for a small elec- 
trode (Fig. [2 top view). As a proxy, we change the num- 
ber of good carbon-metal connections by adjusting the contact 
width, in this way simulating the changing contact quality. 

We consider three widths for the metal electrodes - 3, 5, 
and 7 atomic layers (denoted by 3L, 5L, and 7L) - and two 
kinds of metals - Au and Pd. The electronic states of the two 
electrodes are very different: Au has active s states while Pd 
has only d states. We consider three kinds of CNTs: (5,5), 
(10,0), and Na-doped (10,0) - denoted Na@(10,0) hereafter 
- in which dopant Na atoms are adsorbed periodically on the 
inner surface of the tube. [The adsorption position is deter- 
mined by minimizing the atomic force on the Na atom and is 
found to be slightly off the tube axis and approximately above 
one of the C atoms.] The concentration is one Na atom per 
40 C atoms. Because of the large charge transfer from Na to 
the tube (~0.7 electrons per Na atom by Mulliken population 
analysis), the (10,0) tube is therefore heavily n-doped. 

The electronic structure of the junctions is calculated by 
density functional theory (DFT) using periodic boundary con- 
ditions and a localized basis set with a finite range 18 . In the 
device region, as shown in Fig. ^ the number of atoms in- 
cluded is ~ 350 - 500, depending on the metal width and 
material. The large lateral dimensions (see Fig.^(top view)) 
ensure that the separation among the tube and its images is 
larger than 12A which is much larger than the range of the 
basis functions used. As a result, a good convergence with 
respect to the lateral dimensions can be expected. The con- 



vergence with respect to the length of the tube in the device 
region is also found good, as will be discussed later. 

We adopt a DFT supercell which is larger than the device 
region shown in Fig.^for all the different systems: ~ 500 
- 700 atoms are included, depending on the electrode width 
and material. To avoid too expensive computational cost, a 
single-zeta plus polarization basis set (SZP) is adopted for all 
atomic species. Our test calculation for a small system shows 
that the use of SZP results in only minor differences from re- 
sults using a higher-level double-zeta plus polarization basis 
set. We use optimized Troullier-Martins pseudopotentials 19 
for the atomic cores and the PBE version of the generalized 
gradient approximation 20 for the electron exchange and cor- 
relation. The contact atomic structure is optimized by elimi- 
nating the atomic forces on the atoms around the contact re- 
gion, i.e., the carbon atoms underneath the metal and the metal 
atoms contacting the tube. The relaxation is found quite small 
(see Fig.^(top view) for the (5,5)-Au system) because of the 
choices of the tubes, (5,5) and (10,0), and the orientation of 
the metals, (001). 

Electron transport through the CNT-metal junction is calcu- 
lated using a Green function method^i 2 ^, in which the device 
region and the leads (not shown in Fig.^ are treated exactly 
on the same footing 2 ^. The retarded Green function of the de- 
vice region, Gd(E), is determined by the Hamiltonian given 
by DFT [Hd] combined with the self-energies for the semi- 
infinite leads 
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where Nt is the number of leads or terminals. The transmis- 
sion at any energy, T(E), is calculated from the Green func- 
tion, and the conductance, G, then follows from a Landauer- 
type relation. We adopt the following notation: X^tcim. is the 
two-terminal transmission between tube and metal (terminals 
1 and 3) without electrode 2 attached, which means that the 
tube on the right side of Fig.l (side view) is somehow ter- 
minated and all the electron wave entering that part of tube 
will be totally reflected. T stra ight is the transmission straight 
through the junction from tube to tube, r turn is the tube-to- 
metal transmission with the presence of electrode 2, and R is 
reflection from the tube back to itself. 

In Fig.|2]we show the transmission T2t clm . (E) for the (5,5)- 
Au and (5,5)-Pd junctions with different widths of the metal 



3 




-1 o 1 

E-E F (eV) 



1 o 1 

E-E F (eV) 



FIG. 3: Three-terminal transmission functions T a traight(.E), Tt urn (E), and R(E) for (5,5)-Au [first column, (a)-(c)] and (5,5)-Pd [second 
column, (d)-(f)]. The width of the metal electrode is indicated in the legends. Also plotted is the transmission function of a pure (5,5) tube as 
a comparison (dashed). The Pd traces converge much more quickly than those for Au. 



electrodes. As a comparison, T(E) for the pure (5,5) tube 
is also plotted; it shows, as expected, perfect steps and two 
channels around the Fermi energy. Note that throughout this 
work the position of the Fermi energy is set to be exactly the 
same as in the pure tube. This is required since very far from 
the device region the tube (lead) is completely charge neutral. 

In Fig. |2] the dependence of X^tcrm. (E) on the electrode 
width for the two metals is clear: For Au, the dependence is 
strong. In contrast, for Pd the dependence is much weaker: 
the overall shape of T2tcrm. (E) for the three widths is similar, 
and the result is nearly the same for 5L and 7L. This result for 
Pd also indicates that in our calculation the convergence with 
respect to the length of the tube included in the device re- 
gion is good. Despite the difference in width dependence, the 
equilibrium conductance, G, of both systems increases with 
increasing electrode width: for the thinnest, hence poor qual- 
ity, electrode G ~ Go, while it approaches the physical limit 
of 2 Go for the widest electrode. 

Even the widest electrode (7L) considered here is, of 
course, still much thinner than those used in experiments. 
However, our calculation bears on real experimental situa- 
tions because the number of good carbon-metal connections 
may be similar in that the contact structure considered here 
is perfect while in real experimental situations it is usually 



far from perfect. For example, Au is thought not to wet 
the tube surface but rather will form nanoparticles near the 
surface^. As a result, although the Au electrode can be very 
wide, the tube will pass through the space between these Au 
nanoparticles with very few good C-Au connections formed. 
In this case increasing the diameter of the tube will improve 
the contact quality as more C-Au connections will be formed 
and, therefore, increase the conductance. In fact, a recent ex- 
periment using tubes with different diameters 24 shows that a 
larger diameter yields a larger conductance. Obviously, in- 
creasing the electrode width here is similar to increasing the 
tube diameter: both increase the number of C-Au connec- 
tions and so increase the conductance. A key result here is 
that for good contact quality, the equilibrium conductance ap- 
proaches the physical limit of 2 Go for both Au and Pd. This 
finding differs from a previous first-principles calculation 6 
which used a short segment of tube contacted by two Au elec- 
trodes through several carbon-7r bonds, and also differs from 
a previous model calculation 10 which shows that CNT/metal 
conductance will not be larger than 1 Go, but is consistent 
with recent experiments^. 

The three-terminal transmission in the (5,5)-metal case is 
shown in Fig. [3] for the different metal widths. [Note that re- 
flection R(E)=T tuhe (E)-T stmigH (E)-T tuIn (E).] For the 
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FIG. 4: Two-terminal transmission, T2 te rm.(E), with semiconducting tubes: (a) Na@(10,0)-Au and (b) Na@(10,0)-Pd systems with an 
electrode width of 5L. The transmission function of the bare Na@(10,0) tube is shown in each panel for comparison; in addition, (a) includes 
72term.(£) f° r tne undoped (10,0)-Au system. 



widest case (7L) (see Figs. [3] (a) and (d)) a striking feature 
of the result for Pd is that the straight-through transmission 
is almost turned off while transmission into the metal elec- 
trode dominates. This feature is in very good agreement with 
a recent experiment 11 showing that a Pd electrode suppresses 
inner electron transport through the tube. For Au, T turn is also 
larger than T s t ra ight but they are still comparable and the latter 
is not turned off. 

In Figs. 13 we see that for Au the apportioning of transmis- 
sion among the different components depends significantly 
on the electrode width: For small width, T stra ight (E) and 
Tturn (E) are comparable around the Fermi energy, while 
for the widest electrode (7L), T tuln (E) begins to dominate 
due to the improved contact quality. In the Pd system, the 
convergence as a function of width is much faster (as for 
the two-terminal case Fig. [2j; indeed, one sees only minor 
changes between the curves. This quick convergence indi- 
cates that the carrier injection takes place mainly at the very 
edge of the junction, being consistent with the experimental 
observation 1 '. Thus for Pd, the contact quality is already good 
(i.e., the width is large enough) even though the junction is 
very thin. Since the contact structure is the same, this differ- 
ence between Au and Pd electrodes is due to their different 
electronic states: The Pd d states have stronger interaction 
with C p states than Au s-d states. As a result, Pd is a better 
electrode material than Au. 

Two additional features of Figs. [3] are worth comments. 
First, rt urn (E 1 ) is always smaller than the two-terminal result 
T2term.(-E') of Fig. |2 This is reasonable since in the three- 
terminal case some of the initial flux escapes into electrode 
2 while for two terminals that flux will be totally reflected. 
Second, the reflection increases in all cases near the energies 
where more modes start to propagate in the pure tube. These 
are examples of threshold singularities^. 

We now turn from metallic tubes to discussing the doped 
( 1 0,0) semiconducting tube. The doping is achieved by adding 
Na atoms periodically on the inner surface of the (10,0) tube 
with a high concentration as described earlier. 

The two terminal transmission is shown in Fig. |4] for 
both Au and Pd electrodes. For comparison, the transmis- 
sion is shown for two additional cases: T 2 tcrm.(E) for the 
Na@(10,0) tube alone and for the undoped (10,0)-Au system. 
For the undoped (10,0)-Au system, the Fermi energy is at the 



middle of the gap because one of the leads is the semicon- 
ducting (10,0) tube. After the tube is n-doped, its Fermi en- 
ergy enters the conduction band, and the transmission through 
the Na@(10,0) tube in the absence of any metal has a value 
of 2 around the Fermi energy. As in the metallic case, here 
the Pd electrode results in a larger two-terminal conductance 
(0.79 G ) than the Au electrode (0.56 G ). Both are smaller 
than the two-terminal conductance of the (5,5) systems. 

An interesting feature of the results in Fig. 0] is that 
Taterm^-E 1 ) for the doped Na@(10,0)-Au case is approxi- 
mately a constant shift of that for the undoped (10,0)-Au sys- 
tem. Thus the role of the dopants is mainly to give electrons 
to the tube, and therefore to shift its Fermi energy, without 
noticeably altering the electronic structure of the tube. Our 
calculations indicate this is a general rule of thumb. Follow- 
ing this idea, we can approximately get the two-terminal con- 
ductance of a p-doped (10,0)-metal junction by simply shift- 
ing the Fermi energy slightly into the valence band (to around 
— leV in Fig.|4] For a Pd electrode, this p-type conductance 
is about 1 Go which is in good agreement with a recent exper- 
imental report where the Fermi energy of a semiconducting 
CNT contacted by Pd electrodes was shifted into the valence 
band by applying a back gate [see Fig. 1 (c) in Refill]. 

The three-terminal transmission functions for the semicon- 
ducting tubes cases are shown in Fig. [5] For Na@(10,0)-Au, 
the straight-through and tube-metal transmission around the 
Fermi energy are comparable, just like for the corresponding 
(5,5) system [Fig. 0(b)]. Both of them are, however, much 
smaller than the reflection because the Fermi energy now is 
at the edge of the gap. For a Pd electrode, on the other hand, 
Tturn is comparable to the reflection, near the Fermi energy, 
and both of them are much larger then the straight-through 
transport. Therefore, for an n-doped semiconducting (10,0) 
tube, the Pd electrode turns off conduction through the tube, 
just like for the metallic (5,5) tube. 
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FIG. 5: Three-terminal case with semiconducting tubes: T stra ight (E), T tuI n(E), and R(E) for (a) Na@(10,0)-Au and (b) Na@(10,0)-Pd with 
an electrode width of 5L. For comparison, the transmission function of the Na@(10,0) tube is also shown. 



In summary, we have calculated transmission through both 
2- and 3- terminal CNT-metal junctions using an ideal struc- 
ture in which the tube is surrounded by a Au or Pd electrode. 
In this way we have established the physical limit to which 
experiments will approach by improving the contact quality. 

The main quantities studied are the contact transparency 
in the case of two-terminal systems and the division among 
the different transmission coefficients for three-terminals. The 
two main findings are: 

(1) For two-terminal junctions, when the tube is metallic 
the conductance will approach 2 Go as the contact quality im- 
proves, while for the doped semiconducting tube systems, the 
conductance is about 1 G . In both cases, Pd is better than Au 



for contact transparency. 

(2) For three terminal junctions, the relative magnitude of 
the different types of transmission depends significantly on 
the metal material. The "better" Pd electrode yields near per- 
fect transmission between the tube and electrode and, there- 
fore, turns off the straight transport through the tube for both 
the (5,5) and Na@(10,0) tubes. 

Our results are in good agreement with recent experiments 
and clarify a fundamental discrepancy between theory and ex- 
periment. 
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